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Abstract. Possible phase transition of strongly interacting matter from hadron to a 
quark-gluon plasma state have in the past received considerable interest. It has been sug- 
gested that this problem might be treated by percolation theory. The clustering of color 
sources with percolation (CSPM) is used to determine the equation of state (EOS) and the 
transport coefficient of the Quark-Gluon Plasma (QGP) produced in central A-A collisions 
at RHIC and LHC energies. 
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1 Introduction 

One of the main goal of the study of relativistic heavy ion collisions is to study the deconfined mat- 
ter, known as Quark-Gluon Plasma (QGP), which is expected to form at large densities. It has been 
suggested that the transition from hadronic to QGP state can be treated by the percolation theory (TJ. 
The formulation of percolation problem is concerned with elementary geometrical objects placed on 
a random d-dimensional lattice. The objects have a well defined connectivity radius A, and two ob- 
jects can communicate if the distance between them is less than A. Several objects can form a cluster 
of communication. At certain density of the objects a infinite cluster appears which spans the entire 
system. This is defined by the dimensionless percolation density parameter £ [2|. Percolation theory 
has been applied to several areas ranging from clustering in spin system to the formation of galaxies. 
Figure 1 shows the transition from disconnected to connected system at high densities. 

In nuclear collisions there is indeed, as a function of parton density, a sudden onset of large scale 
color connection. There is a critical density at which the elemental objects form one large cluster, 
loosing their independent existence. Percolation would correspond to the onset of color deconfinement 
and it may be a prerequisite for any subsequent QGP formation. 

The determination of the EOS of hot, strongly interacting matter is one of the main challenges 
of strong interaction physics. Recent lattice QCD (LQCD) calculations for the bulk thermodynamic 
observables, e.g. pressure, energy density, entropy density and for the sound velocity have been re- 
ported [3|. Recently, attention has been focused on the shear viscosity to entropy density ratio rj/s 
as a measure of the fluidity [4 5 6 7 |. The observed temperature averaged 77/5, based on viscous hy- 
drodynamics analyses of RHIC data, are suggestive of a strongly coupled plasma 1 8 9 1 . In this talk a 
percolation model coupled with thermodynamic al relations has been utilized to obtain EOS and the 
rj/s from the experimental data. 



2 String Interactions and Percolation 

Multiparticle production at high energies is currently described in terms of color strings stretched be- 
tween the projectile and target. Hardonizing these strings produce the observed hadrons. The strings 
act as color sources of particles through the creation of qq pairs from the sea. At low energies only 
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Fig. 1. Left: Disconnected system. Right: Connected system 



valence quarks of nucleons form strings that then hadronize. The number of strings grows with the 
energy and with the number of nucleons of participating nuclei. Color strings may be viewed as small 
discs in the transverse space filled with the color field created by colliding partons. Particles are pro- 
duced by the Schwinger mechanisms [ 101 . With growing energy and size of the colliding nuclei the 
number of strings grow and start to overlap to form clusters II 11121 . At a critical density a macro- 
scopic cluster appears that marks the percolation phase transition. 2D percolation is a non-thermal 
second order phase transition, but in CSPM the Schwinger barrier penetration mechanism for particle 
production and the fluctuations in the associated string tension due to the strong string interactions 
make it possible to define a temperature. Consequently the particle spectrum is "born" with a thermal 
distribution ff3l . With an increasing number of strings there is a progression from isolated individ- 
ual strings to clusters and then to a large cluster which suddenly spans the area. In two dimensional 
percolation theory the relevant quantity is the dimensionless percolation density parameter given by 

tnna 

t-ir m 

where N is the number of strings formed in the collisions and S \ is the transverse area of a single 
string and S n is the transverse nuclear overlap area. The critical cluster which spans S appears for 
£ c ^ 1-2 Ifl4l . As £ increases the fraction of S # covered by this spanning cluster increases. 

The percolation theory governs the geometrical pattern of string clustering. It requires some dy- 
namics to describe the interaction of several overlapping strings. It is assumed that a cluster behaves 
as a single string with a higher color field corresponding to the vectorial sum of the color charge of 
each individual string. Knowing the color charge, one can calculate the multiplicity /j.„ and the mean 
transverse momentum squared (pf) of the particles produced by a cluster of strings. One finds 111 11121 



where pi and (pj)\ are the mean multiplicity and average transverse momentum squared of particles 
produced by a single string . In the saturation limit, all the strings overlap into a single cluster that 
approximately occupies the whole interaction area, one gets the following universal scaling law 
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In the limit of high density one obtains 



(Pt >» = 7, J"l ( 4 ) 
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with 



1 -e~S 

being the color suppression factor. It follows that 

» = NF{®p u {p 2 t ) = -±-{p 2 ) l (7) 

A similar scaling is found in the Color Glass Condensate approach (CGO B15I161 . The saturation 
scale Q s in CGC corresponds to (p 2 )i/F(g) in CSPM. The net effect due to F{%) is the reduction in 
hadron multiplicity and increase in the average transverse momentum of particles. The CSPM model 
calculation for hadron multiplicities and momentum spectra was found to be in excellent agreement 
with experiment H17I18I . 



3 Color Suppression Factor F{%) 

The suppression factor is determined by comparing the pp and A+A transverse momentum spectra. 
To evaluate the initial value of F(g) from data for Au+Au collisions, a parameterization of pp events 
at 200 GeV is used to compute the p, distribution II 1912012 II 

dN e /dp* = a/( Po + Pt f (8) 

where a is the normalization factor, po and a are parameters used to fit the data. This parameterization 
also can be used for nucleus-nucleus collisions to take into account the interactions of the strings lfl2l 

dNJdpj = (9) 

(PoyJF(t pp )/F(& + p t ) 

In pp collisions F(g) ~ 1 at these energies due to the low overlap probability. Figure 2 shows a plot 
of F(£) as a function of charged particle multiplicity per unit transverse area ^ /S # for Au+Au 
collisions at 200 GeV for various centralities for the STAR data [19 20 21 1. F(g) decreases in going 
from peripheral to central collisions. The £ value is obtained using Eq. (6), which increases with the 
increase in centrality. These experimental values of F(g) are used to obtain temperature, EOS and rj/s. 



4 Temperature 

The connection between the measured £ and the temperature T(^) involves the Schwinger mechanism 
(SM) for particle production. The Schwinger distribution for massless particles is expressed in terms 
of p] IfTUfol 

dnldp] ~ e-" p - lx2 (10) 

where the average value of the string tension is (x 2 ). The tension of the macroscopic cluster fluctuates 
around its mean value because the chromo-electric field is not constant. The origin of the string fluc- 
tuation is related to the stochastic picture of the QCD vacuum. Since the average value of the color 
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Fig. 2. Color suppression factor F(£) as a function of /S N (fm~ 2 ). The solid red circles are for Au+Au colli- 
sions at 200 GeV(STAR data) 1191 . The error is smaller than the size of the symbol. The line is fit to the STAR 
data. The solid blue squares are for Pb+Pb at 2.76 TeV. 



field strength must vanish, it can not be constant but changes randomly from point to point [13]. Such 
fluctuations lead to a Gaussian distribution of the string tension for the cluster, which gives rise to the 
thermal distribution 113 211 
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dn/dp t ~ e 

with (x 2 ) = n{p 2 )\IF(^). The temperature is expressed as 



(11) 
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5 Equation of State 

Among the most important and fundamental problems in finite-temperature QCD are the calculation 
of the bulk properties of hot QCD matter and characterization of the nature of the QCD phase tran- 
sition. The QGP according to CSPM is born in local thermal equilibrium because the temperature is 
determined at the string level. We use CSPM and thermodynamical relations to calculate energy den- 
sity, entropy density and sound velocity (EOS). After the initial temperature T > T c the CSPM perfect 
fluid may expand according to Bjorken boost invariant ID hydrodynamics 1221 



T dr Cs/ 


(13) 


dT dT ds 


(14) 


dr ds dr 


ds 

— = -Ts/t 
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(15) 


l s 

* = (i + cf)- 


(16) 


dT , 
— s = C 
ds Cs 


(17) 



where s is the energy density, s the entropy density, r the proper time, and C s the sound velocity. 
Above the critical temperature only massless particles are present in CSPM. The initial energy density 
Sj above T c is given by 11221 

dN c I \ 

et = ±-£ (18) 

^ *J pro 
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To evaluate e, we use the charged pion multiplicity dN c /dy at midrapidity and S„ values from STAR 
for 0-10% central Au+Au collisions at yfsim =200 GeV [20|. The factor 3/2 in Eq.(17) accounts for 

the neutral pions. The average transverse mass (m t ) is given by (m t ) = J((pt) 2 + m^), where (p t ) is 
the transverse momentum of pion and mo being the mass of pion. The dynamics of massless particle 
production has been studied in two-dimensional quantum electrodynamics (QED2). QED2 can be 
scaled from electrodynamics to quantum chromodynamics using the ratio of the coupling constants 
Il23l . The production time r pro for a boson (gluon) is IflOll . 



2.405ft 

(m t ) 



T pro = -JZ7T- (19) 



From the measured value of £ and s it is found that e is proportional to £ for the range 1.2 < ^ < 2. 
Figure 3 shows a plot of energy density as a function of e,- = 0.788 £, GeV//m 3 [ 19 20l. 
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Fig. 3. Energy density 6 as a function of the percolation density parameter f . The value for LHC energy is shown 
as blue square. 



Figure 4 shows a plot of s/T 4 as a function of T/T c . The lattice QCD results are from HotQCD 
Collaboration [3]. CSPM is in excellent agreement with the LQCD calculations in the phase transition 
region for T/T c <1.5. The sound velocity and entropy density results are also in good agreement with 
LQCD calculations Iffll . 
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Fig. 4. e/T 4 as a function of T/r c .The lattice QCD calculation is shown as dotted blue line [ 3 1. 
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6 Shear Viscosity 

The relativistic kinetic theory relation for the shear viscosity over entropy density ratio, rj/s is given 
byiH 

H , l^lL (20) 
s 5 

where T is the temperature and A m f p is the mean free path A m f p ~ j^—y n is the number density of an 
ideal gas of quarks and gluons and cr rr the transport cross section for these constituents. In CSPM the 
number density is given by the effective number of sources per unit volume 

AT 

(21) 



S mL 



L is the longitudinal extension of the source, L = 1 fin. The effective no. of sources is given by the total 
area occupied by the strings divided by the effective area of the string S \F{^) IT251 . 

_ (l-e-t)S N 

JV sources n r>sjs\ K^-^J 



j]/ s is obtained from £ and the temperature 

j] TL 



s 5(1 -e-f) 



(23) 



Figure 5 shows a plot of rj/s as a function of T/T c . The lower bound shown in Fig. 5 is given by 
AdS/CFT |27|. The theoretical estimates of rj/s obtained for both the weakly (wQGP) and strongly 
(sQGP) coupled QCD plasma are shown in Fig. 5 JS). It is seen that at the RHIC top energy rj/s is 
close to the sQGP. Even at the LHC energy it follows the trend of the sQGP. By extrapolating the r\j s 
CSPM values to higher temperatures it is clear that rj/s could approach the weak coupling limit near 
'/•//', -5.<S. 
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Fig. 5. r\js as a function of T/J c . Au+Au at 200 GeV for 0-10% centrality is shown as solid black square. wQGP 
and sQGP values are shown as dotted blue and green lines respectively 1 8 1. The estimated value for Pb+Pb at 2.76 
TeV for 0-5% centrality is shown as a solid blue square. The red dotted line represents the extrapolation to higher 
temperatures from the CSPM. The hadron gas value for r]/s ~ 0.7 is shown as solid black circle at T/T c ~0.88 
(261. 
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7 RHIC to LHC 



The STAR results for Au+Au collisions at -\/Jmn = 200 GeV have been extrapolated to estimate F(£) 
values for Pb+Pb collisions at y.sjviv=2.76 TeV. The STAR points in Fig. 2 has the functional form 



Recently, the ALICE experiment at LHC published the charged-particle multiplicity density data as 
a function of centrality in Pb+Pb collisions at y/lim = 2.76 TeV [28 1. The ALICE data points are 
shown in Fig. 2. For central 0-5% in Pb+Pb collisions £ = 10.56 as compared to £ = 2.88 for central 
Au+Au collisions at 200 GeV. The extrapolated value of s for central Pb+Pb collision at 2.76 TeV 
is 8.32 GeV/fm 3 as shown in Fig. 3. For Pb+Pb collisions the temperature is ~ 262.2 MeV for 0-5% 
centrality, which is ~ 35 % higher than the temperature from Au+Au collisions J2TJ. 

One way to verify the validity of extrapolation from RHIC to LHC energy is to compare the energy 
density expressed as s/T 4 with the available lattice QCD results. The LHC point is shown in Fig. 4. 
It is observed that at LHC energy the CSPM results are in excellent agreement with the lattice QCD 
results. The lattice and CSPM results are available for T/T c < 2. 

The estimated value of rj/s for Pb+Pb is shown in Fig. 5 at T/T c = 1.57. These results from STAR 
and ALICE data show that the rj/s value is 2.5 and 3.3 times the KSS bound Il27l . 



8 Summary 

Two central objectives in the experimental characterization of the QGP are the Equation Of State 
(EOS) and the shear viscosity to entropy ratio rj/s. The percolation analysis of the color sources applied 
to STAR data at RHIC provides a compelling argument that the QGP is formed in central Au+Au 
collisions at ysiviv = 200 GeV. It also suggests that the QGP is produced in all soft high energy 
high multiplicity collisions when the string density exceeds the percolation transition. We found rj/s 
= 0.204 ±0.020 at T/T c = 1.15 ( RHIC ) and tj/s =0.260 +0.020 at T/T c = 1.57 (LHC). In the phase 
transition region 77/ s is 2-3 times the conjectured quantum limit for RHIC to LHC energies. The whole 
picture is consistent with the formation of a fluid with a low shear to viscosity ratio. Thus Clustering 
and percolation can provide a conceptual basis for the QCD phase diagram which is more general than 
the symmetry breaking ||29l . 
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